Superconductivity and Antiferromagnetism: Hybridization Impurities in a Two-Band 

Spin-Gapped Electron System 
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We present the exact solution of a one-dimensional model of a spin- gapped correlated electron 
system with hybridization impurities exhibiting both magnetic and mixed-valence properties. The 
host supports superconducting fluctuations, with a spin gap. The localized electrons create a band 
of antiferromagnetic spin excitations inside the gap for concentrations x of the impurities below some 
critical value Xc- When x = Xc the spin gap closes and a ferrimagnetic phase appears. This is the 
first example of an exactly solvable model with coexisting superconducting and antiferromagnetic 
fiuctuations which in addition supports a quantum phase transition to a (compensated) ferrimagnetic 
phase. We discuss the possible relevance of our results for experimental systems, in particular the 
U-based heavy-fermion materials. 

PACS numbers 74.62.Dh, 75.30.Mb, 74.25.Ha, 75.20.Hr 



There has recently been a renewed interest in the role 
of impurities in spin-gapped electron systems. Examples 
include magnetic or mixed- valent impurities in ordinary 
BCS-likc superconductors jl], heavy fermion alloys 
and underdoped high- Tc cuprates as well as disor- 
der driven quantum phase transitions in superconduct- 
ing thin films Q. While many effects from impurities 
in these systems are well understood, such as the pair 
breaking in an ordinary superconductor due to the vi- 
olation of time-reversal symmetry by magnetic impuri- 
ties [||, others remain to be clarified — the coexistence 
of superconductivity and antiferromagnetism in certain 
U-based heavy fermion compounds |^,^ being a prime 
example. 

In this paper we study an exactly solvable model of a 
finite concentration of hybridization impurities in a one- 
dimensional (ID) multiband correlated electron host with 
a spin gap. Our choice of a multiband model is in part in- 
spired by recent experiments on the heavy fermion com- 
pound Ui_2;Th2;Bei3, showing a scattering of electrons 
off impurities in two channels Equally important, 
multiband electron models are intrinsically interesting, 
in particular when allowing for interactions with defects, 
when they generically exhibit non- Fermi liquid behavior 
1^ . Here we address the question how this behavior gets 
influenced by electron correlations in a spin- gapped host. 

The exact solution of our model — to be defined below 
— reveals several remarkable properties. Most interest- 
ingly, for a finite concentration of impurities, too weak 
to destroy the spin gap in the host, there is an impurity- 
induced band of antiferromagnetic spin excitations living 
inside the gap, suggestive of a coexistence of supercon- 
ducting and antiferromagnetic fluctuations. As the im- 
purity concentration increases, the gap closes, signaling 
a quantum phase transition to a fully compensated ferri- 
magnetic phase. To the best of our knowledge, ours is the 
first exactly solvable model where antiferromagnetic fluc- 



tuations of localized electrons and superconducting fluc- 
tuations of itinerant electrons are taken into account si- 
multaneously, allowing for a fully non-perturbative anal- 
ysis of the relevant physics. 

As Hamiltonian for the electron host we take [flOi: 



dxi ~ijj^^{x)d'^llJrna{x) + C dx' 



"Hho 



x5{x - x')iljl^„{x)il;l^,^,{x')il;m'a{x')ipma'{x) 
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where the electron fields ipma carry band indices m= 1, 2 
and spin indices a = ±1/2, with repeated indices in 
(|l|) summed over. The screened electron-electron inter- 
action is here simulated by a local exchange potential 
with strength c > 0, with electrons in the same band 
(and opposite spin) experiencing an attraction while elec- 
trons in different bands (with the same spin) repel each 
other. This interaction, which is integrable |10|, produces 
Cooper-like singlet pairing bound states at any temper- 
ature, but since the model is one-dimensional there is 
no global phase coherence between pairs, and hence no 
long-range order. However, the pairing implies a criti- 
cal field i/c at zero temperature below which there is no 
magnetic response, reminiscent of the Meissner effect (al- 
though there is of course no manifest Meissner effect in 
the system). 

Next, we introduce a Hamiltonian term describing a 
(band-neutral) hybridization of itinerant electrons with 
impurities, with Ni counting their number [|ll|: 



n,mp-ei^J2 J2 [^\SMm){SMm\ + 
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Y,AM,M'{i^lAxj)\SMm){S'Afm\ + H.c.)] . (2) 
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Here \SMm) denotes an impurity state with spin S and 
projection M {S' = S—l/2 > 0) satisfying the constraint 
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Y.ja,\S'M'm){S'M'm\ +Y.M\SMm){SMm\ = 1,6 
measures the ofF-resonance shift of the impurity from the 
Fermi level, ^m_m' = -'2.Sc{Ma\M'f with the Clebsch- 
Gordan coefficients selecting states with M' = M + a as 
{Ma\M') = (S'i;Mcr|5(5 - i);iAf'). Each impurity 
can temporarily absorb the spin of an itinerant electron 
to form an effective spin (S'— i), i.e. the impurity appears 
in two spin configurations jTT|. Note that the spin sec- 
tor of Ti-irnp-ei corrcsponds to a local antiferromagnetic 
exchange 2cSj ■ t{x) between the impurity spin Sj and 
the electron spin density t{x). The two-particle scat- 
tering matrices of the model in (|^) and (||) satisfy the 
Yang-Baxter relations, thus providing for complete in- 
tegrability The model remains integrable for any 

distributions of Sj and 0j but for simplicity we here fo- 
cus on the case where all impurities carry the same spin, 
with an equal number of positive and negative shifts ±9. 

The energies and eigenstates of the model are pa- 
rameterized by three sets of rapidities: charge rapidi- 
ties {kj}jLi (with N the number of itinerant electrons), 
spin rapidities {^a}a=i (with M the number of "down 
spins"), and band rapidities {C/sj/j^i (with n the number 
of electrons in the m = 1 band). (Note that a crystalline 
field will lift the degeneracy in Eq. (Q), with the bands 
becoming unequally populated.) Each eigenstate corre- 
sponds to a solution of the Bethe Ansatz equations, here 
obtained on a periodic interval of length L, 

M n 



N M 

Yl e^f{K ei(A„ - k,) = W e2{Xa - A/) , 

(3) 
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j=i 7=1 

where e„(j/) = {2y inc) / {2y + inc) , j = 1,...,7V, 
a = 1, . . . M, and /3 = 1, . . . , n, and with the energy 



given by E 



To model a RKKY interac- 



U2 

tion in a real system we add an antiferromagnetic cou- 
pling between neighboring impurities While Eqs. (||) 
are insensitive to this addition, it does cause a shift 
of the energy: E E — 27ra;^^' a2s{Xa — 8) where 
X = Ni/L is the impurity concentration, and a„(?/) is 
the Fourier transform of cxp(— n|pc|/2). The structure 
of the impurity-impurity interaction is similar to that of 
Eqs. (^), (^, but with impurity state operators replacing 
electron fields. In the case when all electrons are localized 
it collapses to the well-known Takhtajan-Babujian spin- 
exchange Hamiltonian The fact that energies and 
eigenstates are blind to the spatial distribution of impuri- 
ties is an artifact of integrability. One should note, how- 
ever, that real systems, including certain heavy fermion 
alloys and dopped superconductors 15|, often exhibit 
a large quasi-degeneracy of the states as function of the 



distribution of impurities for stoichiometric compounds. 
Seen from this perspective, the exact degeneracy in our 
model can be turned to an advantage. 

Within the framework of the standard "string hypothe- 
sis" Q , the rapidities satisfying (^) fall into four classes 
in the thermodynamic limit (with continuous distribu- 
tions): (i) real charge rapidities of unbound itinerant 
electrons; (ii) pairs of complex conjugated charge ra- 
pidities representing Cooper-like singlet spin pairs; (iii) 
spin bound states (A-strings); and (iv) inter-band bound 
states (^- strings). Since in ID the most interesting be- 
havior appears in the groundstate (corresponding to the 
low-temperature phases of higher-dimensional analogs) 
we focus on this in the following. 

The groundstate is obtained by filling up the Dirac seas 
of the low- lying excitations, i.e. by populating all possi- 
ble states with negative energies. An analysis shows that 
when a; = only unbound electrons. Cooper-like spin- 
singlet pairs, and inter-band strings of lengths 1 and 2 
have negative energies, and thus make up the Dirac seas 
defining the groundstate. The introduction of a finite 
concentration x ^ of hybridization impurities drasti- 
cally affect the groundstate (and in fact all states that 
are effectively present at a temperature T <C in the 
case that an external magnetic field H is applied): The 
spin levels of the localized electrons now form a low-lying 
band of spin excitations inside the spin gap of the host, 
given by spin strings of length 25*. To see how this comes 
about, we study the integral equations for the dressed en- 
ergies of the excitations (where the "dressing" is due to 
interactions [Q), which, in the case of degenerate bands 
(zero band splitting) can be done analytically. We first 
integrate the equations for the dressed energies of the 
inter-band excitations, as these do not feel the presence 
of the impurities. The resulting set of equations for the 
dressed energies of unbound electrons (e). Cooper- like 
pairs (ip), and spin 2S'-strings (025) has the form: 

{1 - Gi) * e = k'^ - fi - — - Go * ijj + a2s * (f>2S , 

^-2(A2-^-^)-Go*e-x7r^a2s(A±(?) , 



W * (j>2s = SS'iJ - xV — * £ 



(4) 



where fi is the chemical potential, * denotes the convo- 
lution over the appropriate Dirac seas, G„(y) and W{y) 
arc the Fourier transforms of exp(— n|pc|/2)/2cosh(pc/2) 
and coth(pc/2)[l — cxp(— 2S'|pc|)], respectively, and 
V{p) = cos{pO)coth{pc/2)[Go{p)-G4:s{p)]. As revealed 
by (^, for H = and for sufficiently small concentra- 
tions of the impurities, the Dirac sea of spin-2S' strings 
is filled completely. The value of the spin gap (the gap 
for unbound electron excitations) is the smallest energy 
required to depair the Cooper-like spin-singlet state. The 
gap renormalizes in the presence of a small finite concen- 
tration of hybridization impurities, and we find that 
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dX s scch(7rA/c) 



cS 



[Q ± ey + {csy 



-[(2/7r) tan-i sinh(7rg/c) - 1]) , (5) 



where Q is the Fermi level for paired electrons. The spin 
gap decreases due to localized electrons (with their num- 
ber given by (iVi/7r) tan-i[(Q ± e)/cS]) for large 6 
and closes at a critical concentration Xc (i-c. A(a;c) = 0). 
As seen from (||), the smallest value of Xc is obtained 
when S = 1/2, with Xc increasing with increasing S. The 
gap persists for larger impurity concentrations when the 
number of paired electrons (oc Q) is large. However, and 
this is important, there are no additional unbound elec- 
tron excitations appearing when the gap is open (i.e. for 
X < Xc)- Hence, the presence of the hybridization impu- 
rities does not lead to a pair-breaking for these concentra- 
tions, in contrast to the suppression of superconductivity 
in ordinary BCS-like superconductors |5|. The presence 
of the spin 2S-strings in the gapped phase indeed suggests 
a coexistence of antiferromagnetic and superconducting 
fluctuations. It is important to note that the charge 
and magnetic subsystems are effectively disconnected for 
a; < Xc, as revealed by Eqs. (^): The low-energy (con- 
formal) limit corresponds to a direct sum of free bosons 
(hard-core pairs) with scaling dimension rj = I, and an- 
tiferromagnetic spin strings with rj = r{r + 2)/4(S' + 1), 
r — 1,2,... (spinons with the minimal rj = 1/2 when 
5 = 1/2)0. 

The gap, and hence the energies of unbound electrons, 
may become negative when x > Xc- Thus, a Fermi sea 
for unbound electron excitations (with kp (x y/x — Xc) 
opens up in the absence of a magnetic field, signaling a 
quantum phase transition. For H = there are no holes 
in the distribution of spin 2S'-strings, and hence the total 
magnetization of the system remains zero: The appear- 
ance of itinerant electrons correlates with the number of 
localized electrons, with both particle numbers scaling 
with y^x — Xc, thus producing a compensated ferrimag- 
netic phase of the total system. (Note that there are 
two kinds of magnetic excitations in the system: spin 
25-strings and unbound electrons which carry spin 1/2.) 
This compensating effect is due to the RKKY type in- 
teractions between the impurities (which were absent in 
the attractive Hubbard model with hybridization impu- 
rities recently studied in [l3| , where instead a ferromag- 
netic phase appears). However, Cooper- like pairs are still 
present for x > Xc, reminiscent of gapless superconduc- 
tivity, where the gap closes before superconductivity is 
destroyed ||]. We here note that this is different from 
the supersymmetric t—J model with hybridization impu- 
rities |l^], since for that model all low- lying excitations 
are gapless also with no impurities present. 

Let us briefly discuss the model in the presence of a 
magnetic field H . For H smaller than some critical field 
Hc{x), the spin gap persists for x < Xc, and there are 



no unbound electrons in the system. The critical field 
Hc{x) decreases with increasing concentration x and van- 
ishes when X approaches Xc- Thus, as expected, the spin 
gap decreases with the growth of a magnetic field and 
the concentration of impurities. The critical line sep- 
arating the gapless and gapfuU (ferrimagnetic) phases 
manifests itself in the van Hove singularity of the open- 
ing of the band of unbound electron excitations. The 
critical behavior at the quantum phase transition point 
X — Xc, H = also shows up in the scaling dimensions: 
Since an additional band of unbound electron excitations 
now appears, the dressed charge matrix |]l^ gets enlarged 
from 2 X 2 to 3 X 3 with the off-diagonal components pro- 
portional to '^x — Xc- The emergence of 2S'-spin strings 
when S > 1/2, together with the fact that their band- 
width is relatively narrow, oc x, see Eqs. (|), for this 
case imply the appearance of an additional critical mag- 
netic field Hci at which the band of spin strings becomes 
empty (provided that the concentration of impurities is 
finite, since for x — Hd — He). In addition, there is 
the possibility of an additional quantum phase transition 
for H = at Xci > Xc, where the band of spin strings 
may become empty, with the impurities spin-polarized, 
analogous to the case of the attractive Hubbard chain 
with noninteracting hybridization impurities [T^ ]. The 
additional band (of spin strings) increases the effective 
mass of the electrons, producing a larger coefficient for 
the low-temperature specific heat (oc T). In contrast, 
at the critical lines, the van Hove singularities of empty 
bands (of unbound electrons or spin strings) produce a 
Vt behavior of the specific heat. 

To summarize, in this paper we have studied an ex- 
actly solvable model of spin— S* hybridization impurities 
embedded into a two-band correlated electron host with a 
spin gap. A finite fraction of the electrons localize, pro- 
ducing a mixed valence of the impurities. These mag- 
netic impurities interact, giving rise to a band of low- 
lying antiferromagnetic excitations appearing inside the 
spin gap of unbound conduction electrons. The simulta- 
neous presence of Cooper pairs and spin density waves 
in the spin gap suggests that superconducting and anti- 
ferromagnetic fluctuations in fact coexist in this system. 
This is an important result, more so since it is obtained 
by an exact method taking all possible fluctuations into 
account, with no a priori assumption of a (local) sym- 
metry breaking. We predict that the spin gap closes for 
concentrations of impurities above some critical value, 
at which unbound itinerant (and localized) electrons ap- 
pear. However, the total magnetization of the system 
remains zero for H = 0, implying a quantum phase tran- 
sition to a compensated ferrimagnetic phase. A nonzero 
magnetic field less than the critical field He docs not 
destroy the coexistence of superconducting and antifer- 
romagnetic fluctuations (provided that also the impurity 
concentration is smaller than the critical concentration 

Xc). 
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Finally, let us briefly discuss the possible relevance of 
our results for real multiband systems with hybridiza- 
tion impurities. It is clear that the one-dimensionality 
of the model — required for its exact solvability — in- 
troduces features not seen in higher dimensions. In par- 
ticular, spin-charge separation, implying no response of 
the paired system to an external field coupled to the 
current operator, is believed to be an intrinsically one- 
dimensional phenomenon. On the other hand, our model 
shares a number of important characteristics with real 
materials such as the U-based heavy fermion compounds, 
the alloy Ui^ajTh^jBeia [||-^ being a case in point. Com- 
mon properties include (i) an enhancement of the effec- 
tive mass [0; (ii) non- Fermi- liquid behavior in the nor- 
mal (non-gapped) phase jl^]; (iii) power-law behavior of 
the low-temperature specific heat in the presence of a 
spin gap ; (iv) scattering of the electrons off two con- 
figurations of the impurity ion (v) the presence 
of low-lying magnetic excitations in the superconducting 
as well as in the normal state (vi) coexistence of 
superconducting and antiferromagnetic fluctuations for 
some concentrations of magnetic impurities (U ions) [^,0; 
(vii) a (ferrimagnetic) phase with a weak ferromagnetic 
moment for a certain range of impurity concentrations 
pof ; and (viii) (quantum) phase transitions driven by a 
change of the impurity concentration |^,^ 20 1. Consider- 
ing these similarities, it is tempting to speculate that 
our microscopic model supports the phenomenological 
Ginzburg-Landau approach of for the coexistence of 
superconductivity and magnetism in this class of heavy 
fermion alloys Still, the presence of possible nodes 
on the Fermi surface of unbound electrons in 2D or 3D 
may change the picture |23[ | , as may the fact that in our 
model the mechanism for superconductivity lives outside 
the magnetic ions — in contrast to the conventional pic- 
ture of the U-based heavy fermion superconductors. It is 
interesting to note that while integrability imposes cer- 
tain formal restrictions on possible hybridization terms, 
our particular choice in (2) shares some similarities with 
that recently proposed for Ui-^rThjjBeia on purely phe- 
nomenological grounds [^ . This compound shows an 
admixture of two low-energy configurations of U ions: 
the magnetic state Fg of the 5/^ configuration of U'^+ 
and the non-magnetic doublet Fs of the 5/^ configura- 
tion of U^+. This should be compared with the spe- 
cial case of 5 = ^ in (2), with 9 chosen to remove the 
degeneracy between magnetic and non-magnetic states: 
9 = 4[i?(5/'^) — E{5p)], although in our case the mag- 
netic state is four-fold degenerate, in contrast to ||l^ 
where the magnetic state has a two-fold degeneracy. 

Results on the effect of random distributions of impu- 
rity parameters in spin-gapped electron systems, as well 
as a detailed analysis of the present model will be re- 
ported elsewhere. 
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